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The lightest Kaluza-Klein particle (LKP), which appears in the theory of universal extra dimensions, is one of 
the good candidates for cold dark matter. The gamma-ray spectrum from annihilation of LKP dark matter 
shows a characteristic peak structure around the LKP mass. We investigate the detectability of this peak 
structure by considering energy resolution of near-future detectors, and calculate the expected count spectrum 
of the gamma-ray signal. In order to judge whether the count spectrum contains the LKP signal, the x squared 
test is employed. If the signal is not detected, we set some constraints on the boost factor that is an uncertain 
factor dependent on the substructure of the LKP distribution in the galactic halo. Detecting such peak structure 
would be conclusive evidence that dark matter is made of LKP. 


1. Introduction 

At present, most of the matter in the Universe is 
believed to be dark. The existence of non-luminous 
matter, so-called dark matter, was suggested by F. 
Zwicky in 1930s [i|. The dark matter problem is 
one of the most important mysteries in cosmology 
and particle physics 0- One feasible candidate for 
dark matter is the weakly interacting massive particle 
(WIMP), which appears in the theory of beyond stan¬ 
dard model. WIMPs are good candidates for cold dark 
matter (CDM), where cold implies a non-relativistic 
velocity at the decoupling time in the early Universe. 
CDM comprises a large percentage of the matter den¬ 
sity in the Universe [3|, and is necessary to form the 
present structure of the Universe. 

The theory of universal extra dimensions (UED) 
is a popular theory beyond the standard model Q, 
where universal means that all fields of the standard 
model can propagate into extra dimensions. New par¬ 
ticles predicted by this theory are called Kaluza-Klein 
(KK) particles. Here, we consider the theory of UED 
containing only one extra dimension. The extra di¬ 
mension is compactified with radius R. At tree level, 
the KK particle mass is given by Q 



where n is a mode of the KK tower, and toew is a 
zero mode mass of an electroweak particle. 

We assume that the lightest KK particle (LKP) is 
a feasible candidate for dark matter, and we denote 
it Then, jg the hrst KK mode of the hy¬ 

percharge gauge boson. Dark matter should be elec¬ 
trically neutral and stable particles. Hence, LKP ei¬ 
ther does not interact with the standard model parti¬ 
cles or only weakly interacts with them. In addition, 
LKP should have a small decay rate to survive for a 
cosmological time. In the extra dimension, the four¬ 
dimensional remnant of momentum conservation is to 
be conserved as KK parity, and so LKP is stable. This 


hypothesis corresponds to the LKP mass mg(i) being 
in the range 0.5 TeV < m^^i) < 1 TeV using the above 
condition for CDM density Q. In this paper, we as¬ 
sume the TOb(i) is 800 GeV firstly, then we consider 
the change of result in the mass range of 500 GeV to 
1000 GeV. 

There are many LKP annihilation modes which con¬ 
tain gamma-rays as final products. These include 
gamma-ray “lines” from two-body decays, and “con¬ 
tinuum” emission. Branching ratios into these modes 
can be calculated for B^i pair annihilation [1, @ 

and are not dependent on parameters other than 
This paper considers three patterns for the 
continuum: B^^'> pairs annihilate into (i) quark pairs, 
(ii) lepton pairs which cascade or produce gamma- 
rays, or (iii) two leptons and one photon {l'^l~j). The 
gamma-ray spectrum of the continuum component is 
reproduced in FiglT]as per Ref.j^. In this figure, the 
solid line shows the total number of photons per B^^^ 
pair annihilations, the dotted line shows the num¬ 
ber of photons via quark fragmentation, the dashed 
line shows the number of photons via lepton frag¬ 
mentation, and the dot-dashed line shows the num¬ 
ber of photons from the l'^l~'f mode as a function of 

X = . 

When B^^'i pairs annihilate into photon pairs, they 
appear as a “line” at the to^(i) in the gamma-ray spec¬ 
trum. This is the most prominent signal of KK dark 
matter, while in most theories line models are loop- 
suppressed and thus usually subdominant Q- Thus, 
this study focuses on the detectability of this “line” 
structure by near-future detectors accounting for their 
hnite energy resolution. 

The distribution of dark matter is expected to be 
non-uniform in the Universe, and to be concentrated 
in massive astronomical bodies due to gravity. The 
gamma-ray flux from annihilation of dark matter par¬ 
ticles in the galactic halo can be written as ii 

= (Astrophys) x (Particle phys) 

= J(V') X (2) 
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Figure 1: (Color online) Gamma-ray spectra of 
continuum components. The lines show the number of 
photons multiplied by = (E -,as follows: the 
solid line shows the total number of photons per 
annihilation, the dotted line shows the number via quark 
fragmentation, the dashed line shows the number via 
lepton fragmentation, and the dot-dashed line shows the 
number from the Nl~'y component. We have assumed 
= 800 GeV and mass splitting is 5% at the first 
KK level. 


where the astrophysics factor, represented by a dimen¬ 
sionless function J{ip), is calculated as follows: 


JW = 


1 


1 


8.5kpc \0.3 GeV cm 


-3 


X / p^{l)dl{il)) 


(3) 


/ l.O.S 


The function p{l) is the dark matter density along the 
line-of-sight /(?/'), where ip is the angle with respect 
to the galactic center. The particle physics factor is 
written as 


= Const X N-f{av) 


(4) 


where is the number of photons created per anni¬ 
hilation, and {av) is the total averaged thermal cross 
section multiplied by the relative velocity of particles. 
The value of (tru) is accurately computed for a given 
dark matter candidate, so its uncertainty is small in 
terms of considering the cross section containing only 
an s-wave. However, this is not always the case, be¬ 
cause some models have velocity-dependent cross sec¬ 
tions ©[HI- In addition, the astrophysics factor is 
highly dependent on the substructure of the dark mat¬ 
ter distribution in the galactic halo along the line-of- 
sight. Thus, we should consider the so-called “boost 
factor” which indicates the relative concentration of 
the dark matter in astronomical bodies compared with 
some benchmark distributions. The boost factor is af¬ 
fected by {av) and p^{l), and is defined by 


{p\l)) 


AV 


{pm 


AV 


{av)y~. 


^disp 


{av)y:^y 


(5) 


AV 


Bt. 


= BpXBa 


where 'Cdisp is the velocity dispersion, up is the typical 
velocity at freeze-out, a volume AV is a diffusion scale, 
and /Oo(0 ^ typical dark matter density profile, such 

as Navarro-Frenk-White (NFW) [l^l- Bp could be as 
high as 1000 when accounting for the expected effects 
of adiabatic compression [l3l |. 

In the case of gamma-ray flux from LKP annihila¬ 
tion, the particle physics factor is almost fixed for a 
given model, and the boost factor mostly depends on 
the astrophysical contribution. In this paper, we vary 
only Btot as a parameter which describes our limited 
knowledge regarding the astrophysical contribution, 
and consider constraints on its value from observa¬ 
tion. 

Recent progress in gamma-ray observation has re¬ 
vealed new findings in the galactic center region. The 
high energy gamma-rays from the galactic center have 
been observed by the High Energy Stereoscopic Sys¬ 
tem (HESS) [3, the Large Area Telescope on board 
the Fermi Gamma-Ray Space Telescope (Fermi-LAT) 
[T^ and other experiments. However, the observed 
gamma-ray spectrum is represented as a power-law 
plus an exponential cut-off, and is hardly compatible 
with a dark matter signal [l^. Recently, some ev¬ 
idence regarding the enhancement of the continuum 
component of the gamma-ray emission from the galac¬ 
tic center region detected by the Fermi-LAT has been 
reported and argued as a possible signal from the de¬ 
cay of dark matter particles [Tj, [l^. Further, an 
enhancement of around 130 GeV in the energy spec¬ 
trum of gamma-rays from the galactic center region 
has been reported which may indicate a possible dark 
matter signal [1^ [2l|, [2^ . Discussion relating to 

unifyin g th e continuum and the line has also been pre¬ 
sented [^. However, the analysis by the Fermi-LAT 
collaboration did not confirm the significance of the 
line detection [2^ [2^. Thus, the situation is still un¬ 
clear and more sensitive observation is necessary to 
resolve the issue. 

In the following, we focus on gamma-ray observa¬ 
tion with near-future missions, such as the Calorimet¬ 
ric Electron Telescope (CALET) CALET is a fine 
resolution calorimeter for cosmic-ray observation to be 
installed on the International Space Station. CALET 
will detect gamma-rays in the energy range of 4 GeV 
to 1 TeV with about 1000 cm^ effective area, and a 
few percent enerey resolution, suitable for gamma-ray 
line detection [^ . 

In this paper, we analyze the gamma-ray spectra 
from bO') pair annihilation accounting for the finite 
energy resolution of gamma-ray detectors and pur¬ 
posefully discuss the observability of the “line” at the 
m^(i). We then give possible constraints on the boost 
factor by near-future detectors. 
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2. The effect of energy resolution 


The gamma-ray spectrum d^{Afl)/dE-y reaching 
a detector can be expressed as Q 

~ Const X Btot X (6) 

doc 

where Afi is the angular acceptance of the detector, 

{av) 


Const ~ 3.5 X 10 ® 


3 X 10“2®cm3s“^ 


V rngii) J 


An, (7) 


and (Jgc)ao ^ dimensionless line-of-sight integral 
averaged over An. If we assume an NFW prohle, 
{Jgc) equals to 0.13 for a An = 10“® [2^. In 
this case dNj/dx includes both the continuum and 
line components. 

Now, we discuss the effect of energy resolution. If 
the measured energies of detected gamma-rays behave 
like a Gaussian distribution and the energy resolution 
is 1%, the measured “continuum” gamma-ray spec¬ 
trum is blurred and should appear as shown in Fig|^ 
Here we draw the curve assuming the following equa¬ 
tion 


g{E) oc 


fiE') X exp 


{E-E')^' 


dE', 


( 8 ) 


where f{E') corresponds to a function shown by the 
solid line in FiglTJ and as is the energy resolution. 



Figure 2: (Color online) Gamma-ray spectra of the 
continuum accounting for energy resolution assuming 
= 800 GeV. The solid line assumes an energy 
resolution of 1% with a Gaussian distribution, and the 
dotted line does not include the effect of energy 
resolution, as per the solid line in Fig[T] The assumed 
boost factor is 1000. 


Next, we analyze how the “line” from the pair 
annihilation into photon pairs looks above the “con¬ 
tinuum”. The three patterned lines shown in Fig 13] 



Figure 3: (Golor online) Gamma-ray spectra of 
continuum plus line diffused by the energy resolution 
assuming mg(i) = 800 GeV. The solid line shows the 
continuum component only, assuming the energy 
resolution of 1%, while the dotted, dashed and 
dot-dashed lines show the continuum plus line 
components assuming energy resolution values of 0.5%, 
1%, 2% respectively. The assumed boost factor is 1000. 


assume different energy resolutions which we take as 
0.5%, 1% and 2% with the Gaussian distribution. In 
Fig Is] the solid line shows the continuum component 
only with an energy resolution of 1%, and the pat¬ 
terned lines show “line” plus “continuum” spectra for 
different energy resolutions: the dotted line, dashed 
line and dot-dashed line show the spectra when the 
energy resolution is 0.5%, 1%, 2% respectively, as¬ 
suming the boost factor Btot = 1000. 

We can transform the spectra into counts to be ob¬ 
served by gamma-ray detectors. This is accomplished 
through multiplying by a factor of 0.03 for an assumed 
observation time of 1 yr and an assumed effective 
area of 1000 cm^. These values arise from the typ¬ 
ical aforementioned CALET sensitivity [2^. When 
analyzing observational data, the energy bin width 
must be specified. A bin width of 1% of mg(i) (about 
one standard deviation of energy reconstruction) was 
used in order to avoid energy information loss. The 
resulting histograms are shown in Figjd] where plots 
of the three cases corresponding to energy resolutions 
of 0.5%, 1% and 2% are shown. The figure shows that 
if the energy resolution of the detectors becomes 2% 
or worse, the characteristic peak indicating the 
will be diffused, making it hard to resolve into the line 
and continuum components. Thus, the energy resolu¬ 
tion for gamma-ray detectors should be better than 
2%, in order to “resolve the line” without the need for 
detailed analysis. 

Thus far, we have taken the LKP mass to be 
ixibG) = 800 GeV, and calculated count spectrum 
for its mass. Now, we vary the mass from 500 GeV 
to 1000 GeV in 100 GeV intervals, and calculate the 
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Figure 4: (Color online) Expected count spectra near the 
peak assuming energy resolutions of 0.5%, 1% and 2% 
assuming mg(i) = 800 GeV. The bin width of histograms 
is 8 GeV, equaling 1% of the m^(i). The assumed boost 
factor is 1000. 


count spectrum for each mass. The results are shown 
in Figini This figure shows that the characteristic 
peak structure is visually clearer when to^ci) is heav¬ 
ier. That is, the line component becomes relatively 
larger since the continuum component decreases for 
heavier mg(i). 



Figure 5: (Color online) Expected count spectra, 
assuming energy resolutions of 0.5% and 1%. The data 
space is 8 GeV, which is 1% of for the mg(i) = 800 GeV. 
The assumed boost factor is 1000. 


3. Discussion 

We now discuss the observability of the LKP sig¬ 
nal in near-future detectors, taking account of the ob¬ 
served background spectrum. That is, we give esti¬ 
mates for the accessible range of the boost factor when 
the observed counts are significantly different from the 


background spectrum. Here, we assume the gamma- 
ray spectrum from HESS J1745-290 located near the 
center of the Galaxy is the source of the background. 
Its spectrum is given by 0 

cn \ -2.10±0.04±0.10 

4 ) 

E 

'' [ (15.7±3.4±2.5)TeV_ 

xlO"® TeV^ m-2 s"^ (9) 

Note that with the energy resolution of HESS (15%), 
the LKP “line” signal is broadened and hard to detect. 

To investigate the detectability quantitatively, we 
employ a y-squared test method to judge whether the 
excess counts are statistically meaningful. 

First, we define as 


d<^ 

— = (2.55 ±0.06 ±0.40) 
uE 


X 


E 


([count ± background] — background)^ 
background 


( 10 ) 


where N is the number of energy bins, corresponding 
to degrees of freedom for the y-squared test. We then 
specify the energy range: 


Energy range = [100 GeV, 1 TeV] (11) 

with bin width of 0.8 GeV (= 0.1% for = 800 

GeV). Thus, N is about 1000 in this case. The upper 
bound of the energy range under analysis is fixed as 
Ws(i) ± Zue to allow finite energy resolution. Hence, 
at this energy, the degree of freedom is one {N = 1). 
Then, we vary the lower bound of the energy range to 
lower energies. Thus, N gradually increases as we ex¬ 
pand the energy range to lower energies. Eor example, 
N at the peak for 1% energy resolution is 


N [Epeak, rngd) ± Sctb] = 40. (12) 

We investigate the value of boost factor when y^ is 
bigger than some critical value for each N. The re¬ 
lation between N and the upper bound of the boost 
factor is shown in FiglBl where the “peak” on each line 
corresponds to the value when N equals to Eq. (HH). 
Then, Str^ are the energy width limits within 3a from 
the peak. Thus, they are given as 

N at 3a^ = N [Epeak ± 3aE,mBd) + Sctb] (13) 

One can see from this figure that the limit on the boost 
factor does not change rapidly when we include energy 
bins well below the peak. An accessible boost factor 
would be smaller than 500 when N is in the range 30 
- 200. These values of N correspond to being near the 
peak energy for annihilation of LKP. 

We applied similar analyses for other LKP masses. 
The results shown in Pig (7] indicate that the constraint 
for the boost factor is tighter for lighter . In 
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Number of degrees of freedom, N 

Figure 6: (Color online) Expected limits on the boost 
factor for the mg(i) = 800 GeV as a function of the 
number of degrees of freedom of the observed energy 
range. The horizontal solid line shows a boost factor of 
500. The dotted, dashed and dot-dashed lines show the 
boost factor when values are bigger than critical 
values for energy resolution 0.5%, 1.0%, and 2% 
respectively. 



Figure 7: (Color online) Expected limits on the boost 
factor for each mass, assuming an energy resolution of 
1%, as a function of the number of degrees of freedom of 
the observed energy range. 


addition, we compared the results when using energy 
resolution of 1% and 0.5%, as shown in FigUl The 
number of events near the peak increases with better 
energy resolution, and the resulting constraint near 
the peak is tighter. 

Some constraints from observations on the KK dark 
matter models have been reported. The Fermi-LAT 
team searched for gamma-ray emission from dwarf 
spheroidal galaxies around the Milky Way galaxy 
and set constraints on dark matter models with non¬ 
detection results [S^ . The HESS array of imaging air 
Cherenkov telescopes observed the Sagittarius dwarf 
spheroidal galaxy in the sub-TeV energy region and 



Figure 8: (Color online) Comparison of the expected 
limits of boost factors with 1% and 0.5% energy 
resolution. The individual lines refer to the scenarios 
with masses set from 500 GeV to 1000 GeV in 100 GeV 
intervals. 


derived lower limits on the of 500 GeV [3l|. 

These results put constraints on (av) of dark matter 
halo KK particles. The present limits allow the max¬ 
imum value of boost factors of several to 1.5 x 10'^ 
depending on . Our analysis on future high en¬ 

ergy resolution observation improves the limits on the 
boost factor or the chance to detect the signal. If such 
signals are detected, we will be able to say that dark 
matter is made of LKP, which will be evidence of the 
existence of extra dimensions. 


4. Conclusion 

Energy resolution plays a key role in detecting the 
line structure of the gamma-ray spectrum expected 
from annihilation of LKP dark matter as predicted 
by UED theories. This paper investigated the effects 
of energy resolution of gamma-ray detectors and cal¬ 
culated the expected count spectrum. The predicted 
gamma-ray spectrum is the sum of the continuum and 
a line corresponding to m^ci), but this characteristic 
spectrum is diluted when we account for the finite 
energy resolution of detectors as shown in Figl5] and 
FigO Further, if we assume the exposure (area mul¬ 
tiplied by observation time) of near-future detectors, 
count statistics will be the final limiting factor. The 
characteristic peak indicating the rn^(i) would be dif¬ 
fused if the energy resolution is 2% or worse. How¬ 
ever, with qualitative statistical analysis, we may be 
able to detect a peak statistically by subtracting a 
background from the observed spectrum. In addition, 
if mg(i) is heavy, the observed gamma-ray spectrum 
will show the characteristic peak clearly because the 
continuum component decreases relative to the line 
component. 
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This paper also estimated the accessible range of 
the boost factor using a x-squared test assuming the 
HESS J1745-290 spectrum as a background. If the ob¬ 
served energy range for gamma-rays extends to lower 
energies, the accessible range of the boost factor will 
be lowered since a higher amount of continuum events 
will be detected. If the signal is not detected, the up¬ 
per limit of the boost factor is about 500 if only taking 
data near the peak, and about 100 if the whole energy 
range is covered. Furthermore, if rng(i) is light or the 
energy resolution of the detector is good (say the order 
of 0.5%), we may tightly constrain the boost factor. 

If the gamma-ray line structure is observed in the 
fnture, we may identify LKP dark matter, which will 
provide strong evidence for the existence of extra di¬ 
mensions. 
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